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Abstract 

In this paper we present a number of four-dimensional neutral sig- 
nature exact solutions for which all of the polynomial scalar curvature 
invariants vanish (VSI spaces) or are all constant (CSI spaces), which are 
of relevence in current theoretical physics. 

1 Introduction 

In [1] we considered pseudo-Riemannian spaces of arbitrary signature for which 
all of their polynomial curvature invariants vanish (VSI spaces). We discussed 
an algebraic classification of pseudo-Riemannian spaces in terms of the boost 
weight decomposition and defined the Si- and N-propertics, and showed that 
if the curvature tensors of the space possess the N-property then it is a VSI 
space. We also showed that the VSI spaces constructed possess a geodesic, 
expansion-free, shear-free, and twist-free null-congruence and hence are pseudo- 
Riemannian Kundt metrics of arbitrary signature. 

We will use this fact to construct a set of four-dimensional (4D) neutral 
signature metrics, in which all of the scalar curvature invariants vanish, gen- 
eralizing the 4D degenerate Kundt metrics [3] and the VSI metrics [2] in the 
Lorentzian case. We will also consider 4D neutral Einstein metrics which for 
which all of the curvature invariants are constant [SI [B] (CSI spaces), and present 
two simple examples for illustration. 

All of these 4D neutral signature solutions are of interest in the twistor 
approach to string theory [7] and for spaces admitting parallel spinors [5] (where 
some VSI spaces are already known). 
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2 Pseudo-Riemannian Kundt metrics 

Let us first define the pseudo-Ricmannian Kundt metrics in arbitrary signature 
and arbitrary dimension [T]. Consider a null- vector £ with an associated null- 
vector n so that i^n^ = 1. We then form the null- frame {i,n,u> 1 } and define 
L l i = ^u'w/i where u)^ is the vielbein corresponding to lj 1 and, as usual, 
a semi-colon denotes covariant differentiation. We call £ geodesic if and only if 
(iff) l^i v -^ = 0, twist-free iff Luji = 0, expansion-free iff L t i = 0, and shear-free 
if %) = 0. 

Definition 2.1. A pseudo-Ricmannian space is a Kundt metric if it possesses 
a non-zero null vector £ which is geodesic, expansion-free, twist-free and shear- 
free. 

We will consequently consider such metrics in null coordinates: 

ds 2 = 2du [dv + H(v, u, x c )du + Wa(v, u, x c )dx A ] + qab (w, x c )dx A dx B ', (1) 

where the range of A is over n — 2 coordinates The metric (HJ possesses a null 
vector field £ obeying 

^-,u = Ln£p£ v + Lui^Lu 1 ^, 

and consequently, is geodesic, non-expanding, shear-free and non-twisting. 
In particular, 

= ^ = ^k^) = ^'"^M = °- ( 2 ) 
2.1 Walker metrics 

A special case of the Kundt metrics is the Walker metrics. If £ x and £ 2 are two 
mutually orthogonal null vectors spanning a 2-dimcnsional invariant null plane, 
then they satisfy the recurrence condition 

[*\a%]];c = (3) 

for some recurrence vector k c (and the two null vectors are automatically sur- 
face forming). We can interpret this as requiring that the bivector, £ l A £ 2 
of the invariant 2-plane is recurrent. A pseudo-Riemannian space admitting a 
2-dimcnsional (or 1-dimcnsional) invariant null plane is called a Walker space 

a- 

Furthermore, Walker [J] showed that the metric can be written in the canon- 
ical form: 

d.s 2 = du I {28 IJ dv J + B u du J + H H dx l ) + ,1, .dr'dH. (4) 

where Bjj is a symmetric matrix which may depend on all of the coordinates, 
while Hn and Aij do not depend on the coordinates v 1 . This Walker metric 
admits a null-vector £ x (in fact two, £ x and £ 2 ) which is geodesic, expansion- free, 
shear-free and twist-free, and is hence a pseudo-Riemannian Kundt metric. 
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3 4D Neutral-signature case 

A space is VSI if it satisfies the Si-property and the N-property [T], and as- 
suming a Kundt metric (JlJ , this implies that the Riemann components i?oioA — 
— (1/2)Wa,w = and -Roioi = cr = 0, and the matrices 



R AB 



CD 



in eqns. (20)- (22) for the components of the Riemann tensor in [T] are nilpotent 
(i.e., they have only zero-eigenvalues). In particular, assuming Wa,w = 0, these 
quantities are related to the Riemann components as follows: Rqiab = ^ab, 
Rq(a\i\b) = (1/2)sab and R OD = R CD . 

Let us consider the 4D neutral (2, 2)-signature case. Here we can write 

ds 2 = 2 (I 1 ® n 1 + I 2 (g> n 2 ) . (5) 

There is only one independent component of the "transverse Riemann tensor" 
R AB CP and so, requiring the N-property, implies that this must be flat space. 
Therefore, we can write: 

2£ 2 ® n 2 = 2du 2 <g> dv 2 = -dT ® dT + dX <g> dX. 

We can then find two classes of pscudo-Ricmannian Kundt VSI metrics, 
which can be written: 

ds 2 = 2dw 1 <g) (dv 1 + Hdu 1 + W^dx^ 1 ) + 2du 2 <g> dv 2 , (6) 

where the 1-form W^ v 1 can be null, or timelike/spacelike. The two classes are: 

Null case: 

W^dx" 1 = v^!^ (u 1 ,u 2 )du 2 + W^\u\u 2 ,v 2 )du 2 + w£\u\ u 2 , v 2 )dv 2 , 
H = v 1 H^(u\u 2 ,v 2 ) + H^(u 1 ,u 2 ,v 2 ), (7) 

Spacelike/timelike case: 

W^dx" 1 = v l W^6X + W^ 0) (u\T,X)dT+W { ° ) (u 1 ,T,X)dX, 
(v 1 ) 2 < -as 



H = 
and 



(V (1) ) +v 1 H^(u\T,X) + H^(u\T,X), (8) 



W (1) = ~. where e = 0, 1. (9) 
X 

We note that these possess an invariant null-line if = 0, and a 2- 

dimensional invariant null-plane (Walker space) if wffl = for the null case 
(in order for the spacelike/timelike case to possess an invariant null 2-plane, it 
needs to be a special case of the null case.) 
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3.1 Coordinate transformations 

The form of the metric ((SJ) is invariant under the following transformations in 
the null case: 

(u'VW 2 ) = yy+hyy^^yy) (io) 









= H<®-hHW-h v i 




= w$> 




= w$ - hW$ - h u2 




= W$> - h V 2 




,i/ 2 ) = (g^y/g,^ 




= (ff (1) +v)/(v) 2 




= H^/{g, u if 


w$ 






= W<$/g, u > 




= W$>/g tU > 



H 'W = H (o) _ hH (i) _ h 



(hW^A /8 



W$) = - hWM - h x 

w p = wP~h T 



H'W = {g^H^+g„ u ,)/{g^f 
H'^=H^/(g^f 

= wP/g tU , 
W W=wW/g, u i 



(11) 



Similarily, for the spacelike/timelike case: 

(u n ,v n ,X',T') = y,v 1 + h{u 1 ,X,T),u 2 ,v 2 ) (12) 



v'\X',T') = {g^y/g^XiT) (13) 
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4 Exact VSI solutions 

Let us present the exact VSI solutions in the null and non-null cases. 
4.1 Null case 

Consider first the null case using the form of the metric eq. ([7]). Computing the 
Ricci tensor, the conditions for the VSI space to be a vacuum, R a b = 0, yields 
the following system of partial differential equations: 

av £ au A ov z ov z 

+ H (i) J-W® + H^—W [ t> + WVW { ? —W^ 
du 2 v dv 2 u u v dv 2 u 

+ 2W ( ? — W ( ? + — If 1 ! 1 - w { ? —w { ? 
" 2 du 2 " 2 dvldu 2 v2 " 2 du 1 u2 

dv}dv 2 u 2 \ u2 v2 ) u2 1,2 du 2 " 2 
JLw(°> iw(?) _ I ( I { iw(») 



+ —w { 7— WV - - —W K 7 - - — (14) 
dv 2 u2 du 2 v2 2 \dv 2 u2 ) 2\du 2 v ' { ' 



AjW - + i *!_w<? - 1 92 (15) 

<9« 2 2 " 2 dv 2 v 2d{v 2 ) 2 u2 2dv 2 du 2 v2 v ' 



29U 1 11 <9u 2 2 " 2 V u2 J ^ 2 u2 dv 2 u 

2 " 2 du 2 1,2 2 dv 2 du 2 u2 2 d{u 2 ) 2 v2 1 ' 



= A. w w -(w { V) 2 (17) 

The form of the metric ^ is preserved under the transformations (JTUJ), (fTTTl . 
Applying the transformation (JTUJ), without loss of generality we can set = 
0. 

One solution of eqn. (|T7)) is W}i — 0; the resulting solutions are then [§]: 
H^Xu 1 ,u 2 ,v 2 ) = a(u\u 2 )+ /3(u\v 2 ), 

W {0 2 {u 1 ,u 2 ,v 2 ) = 2a(u 1 ,u 2 )v 2 -2 f 0(u 1 ,v i )dv i + 7 (u\u 2 ), 



H^(u\u 2 ,v 2 ) = 2/3(u\v 2 )v 2 / adu 2 -2^(u 1 ,w 2 )u 2 / f3{u\ v 2 )dv 2 + p(u, y) / -/(u,u 2 )du 



;2 



I ^ a {u\u 2 )du 2 -u 2 [ ^I3(u\v 2 )dv 2 +F 1 (u\u 2 )+F 2 {u\v 2 ) 
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In particular, we obtain the special neutral 4D, Ricci-flat, VSI solution of [T] in 
which all of the W's are zero, H 1 = H 2 = 0, and H^[v}, u 2 , v 2 ) = Fi(u 1 , v 2 ) + 
F2(u 1 , u 2 ) in terms of the arbitrary functions F\ and F2. 

Assuming that W 2 7^ 0, eqn. (fTT)) can then be solved as follows: 

W { V = -2 = -2f 

Next we integrate eqn. (|15p with respect to v 2 , which gives 

cm 2 

where /3 is an arbitrary function. We then obtain 

d 



ou z 



This has the solution 



H^=e- J Pe I du 2 + 1 (u\v 2 )^j , 

where I — J fdu 2 = ln(/ _1 ) and 7 is an arbitrary function (and any integrating 
factor that is a function of (u 1 ^ 2 ) can be absorbed into 7). For notational 
simplicity we redefine the arbitrary function (3 as ^(u 1 ^ 2 ) = f J (5/ fdu 2 . We 
can now write down the solution for H' 1 ' in the general form: 

H^(u\u 2 ,v 2 ) = fi{u\u 2 ) + fj{u\v 2 ). 

We can now use eqn. to obtain: 

W$ =S(u\u 2 )v 2 -2f(u\u 2 ) J 1 {u 1 ,v 2 )dv 2 + V (u\u 2 ) 

Equation (fTB"|) now contains only the terms Wji , and W 3 , and hence 

simply puts constraints on the arbitrary (integration) functions. After some 
manipulation we obtain 

(3(u\u 2 ) = l -8(u\u 2 ) - ^p-f( u \u 2 ) + /V(u\«W>«>* 2 + 
2 au 1 J 

Using the remaining coordinate freedom (jllj) we can now choose e(u 1 ) = 0, and 
hence we obtain the solutions 

W^> =5y- 2/ J 7 d« 2 + r, 
W$> = -2/ 
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Substituting these expressions into (fT^|) we obtain the solution 

H {0) = f(J jdv 2 J Sdu 2 -fa J (7) 2 d« 2 ^) (J f du 



ov 
dv 2 



dt ' 2 ) / ??du2 

7 dv 2 ] [ / / fSdu 2 du 2 



d, 2 ^ - lv 2 ( I [ / <5dw' 



f \du 

du 1 J 2 V dii 1 
1 

7 



u 2 



/If / /Mu 2 ) Mu 2 



+ F 2 {u\u 2 ) 



(18) 



4.2 Spacelike/timelike case 

Setting R a b = and using the coordinate freedom (|T2"j) to set 
obtain the following equations: 



2x2 §f H(1) + m ( x2 ik w x ) = (19) 

2* 2 |!^ (1) = ^ ( 2 ^ 2 ^^ W - 2 ^ ) ( 2 °) 
X 2 ^^ = - (2X 2 A^d) _ 2 ^)j (21) 



W j» ) ( 2X 2 A ff (D _ 2W(°A + 2X 2 wf } Affd) + 2X 2 H< 1 ) A W W 

-^(^<») 2 + «^<> ,22, 



8 



A. Alcolado, A. MacDougall, A. Coley and S. Hervik 



Equations (fT9"|) - (f2~T]) do not contain flf"). We can solve those equations first 
for and and then use the fourth equation to obtain H^ '. Note also 

that the first three equations are not independent. Taking the derivative of (|19[) 
with respect to T and using eqn. (|2"Tj) gives eqn. (|2T)j) . 

Integrating eqn. ([T9]) with respect to T gives 

where a(u, A) is an arbitrary function. Replacing TjW in eqn. (|2ip with this 
form gives an equation for Wj^ in terms of an arbitrary function. After applying 
the general identity x 2 f + 2x-^f = x-^j (xf) for any function f(x, . . .), we 
obtain the more compact equation for : 



A change of variables A = X + T , B = X — T is then used. The operator on 
the left turns into a single second order mixed derivative, so that 

^^(xW^)=2X^-a(u,X), (25) 



and thus after integration 

W x = X~{JJ ^^<u,X)dAdB + p{u,A)+ 1 {u,BYj, (26) 

where /? and 7 are arbitrary functions. Since the integrand in this expression 
is not a function of T, integration with respect to A and B can be replaced by 
integration with respect to X using dA = dB = 2dX: 

= 2X j adX -4 J J adXdX (27) 

where integration by parts has been used. It is convenient to redefine our 
arbitrary function as a(u, X) = 2 JJ adXdX and rescale f3 and 7, where we 
keep the same name for notational simplicity. Hence, 

W4 0) = I (£(«> * + T) + 7(«, X - T) - a(u, X)) + || («, A). (28) 

It is straightforward to write down the solution for by using eqn. 
Keeping in mind the redefinition of a, we obtain: 



H 



a(u, A) - /3(u, X + T) - 7 (u, X - T) 



l|(«,I+T)-I^( U ,I-T)]. (29) 
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We now have solved eqns. (fTi?|) and (|2ip . Eqn. (|20p is automatically satisfied 
since it is not independent of the other two, as was discussed earlier. All that 
remains is to solve eqn. (|22|). 

The terms involving derivatives of with respect to X can be simplified 
by using differential identities, by using the previously solved equations and by 
using 

where we have used eqn. (|21[) and the chain rule. Furthermore, 

2X 2 W i " ) ^-H^ + 2X 2 H^^-W^ = 2X 2 ^-(W^H^) (31) 
x dX dX x dX \ x J y ' 

Therefore, we can write down the more concise equation for H^ ': 



f d 2 


d 2 \ 


(»:>) 


(dX 2 


dT 2 ) 





11 S! : <'4W , '" 1| )-^K , )1 <*o 



X \dudX dX v J AdT 2 

For given functions and this equation can be integrated to solve for 



5 Some exact CSI solutions 

We can also consider 4D neutral signature spaces which arc CSI, for which 
all of the curvature invariants are constant [5]. By solving the appropriate 
components of the Riemann tensor equal to constants (see eqns. (15) - (22) 
in pQ), and using similar techniques to those in [5], we can find examples of 
4D neutral CSI metrics. In general the Einstein equations that need to be 
solved are quite complicated. Therefore, let us consider two simple examples 
for illustration. In particular, the neutral signature "Siklos metrics" (in which 
the only non-zero independent invariant is the Ricci scalar, which is constant) 
give rise to relatively simple equations that can be solved completely. 

Lorcntzian CSI spacetimcs are known to be solutions of supcrgravity theory 
when supported by appropriate bosonic fields 0, and it is likely that neutral 
signature CSI spaces are also of physical interest [3 \8\. 

5.1 CSI example 1: 

Let us write the CSI metric as: 



ds 2 = 21 1 ® n 1 - vtp{2KX)dT ® dT + dX <g> dX, 
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where 

i 1 = du (33) 

n 1 = dv+ [v 2 a + vH {1) (u,T,X) + H (Q} (u,T,X)}du (34) 

+ [Av + W [ x ] (u, T, X)]dX + [Bv + WP (u, T, X)} exp(AA)dT 

There are many cases which lead to an Einstein space. In particular, we 
consider the simple case: 

A = -2K, B = 0, a = 0. 

Writing down the Einstein conditions, we obtain A = — 3K 2 from the diagonal 
terms of the Ricci tensor, and differential equations from the off-diagonal terms. 
All of these have simple curvature structure; in particular, C a b c dC abcd = 0. 
Explicitly, the metric has the form 

ds 2 = 2du ® (dv + [v 2 + vH {1) + H m ]du 

+ [-2Kv + W { x ] ]dX 

+ W4 0) exp(AA)dT) - exp(2KX)dT <g> dT + dX ® dX (35) 

where H arc functions of u,T,X. The conditions for the 

metric to be an Einstein space yield the following system of partial differential 
equations: 

4-^*4 } cxp(2AA) + 2Acxp(2AA)^ r 
+ 4A 2 exp(2AA)ff(°> + 2^- 
° T exp(AA) x 



dX iy ' dT 

>wl 

du 



2if exp(2AA)^-^- - (w| 0) ) 2 ^ 2 cxp(2XA) 



6Kexp(2KX)v^-^- + 2Acxp(2AA)£T (1) w4 0) 

(0) 



2W^KeMKX)^- 

2i/«^ cxp(2AA) - 4^wf exp(AA) 

2H^^eMKX) 
d 2 W<t 0) .„ vs n d 2 H^ 

d 2 H^ 




dx2 x - eM2KX) = (36) 
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BWf d 2 W^ , , d 2 wf , 
K ^t-^*MKX) + 1 ^§ F =* (37) 



2^cxp(2A*)-^f e M KX) 
— K ™ 



exp(KX) + -^f- = (38) 



dT 

The v dependence in eqn. (|3f)|) gives 



Q2 F (1) q H (1) 

- 2 ^p" exp(2KX) - 6Kcxp(2KX)-^- 
d 2 H^ 

+ 2 ^^ = ° (39) 

The form of the metric Q35p is preserved under the transformations ()12[) . 
([TU)) . Thus, without loss of generality, we can set = 0. We note that 

does not appear in equations (|3"Tj) . (|55|) . (|3^|) . We first use eqns. (|3"T)) and ([551) 
to solve for W-jp and H^ l \ and then use eqn. (J39J) to put constraints on these 
solutions. 

First we can integrate eqn. ((57)) with respect to T to get 

2tf« + fcl4 0) + ^-+ a = 0, (40) 
oX 

where a = a(u,X) is an arbitrary function. Now, using using eqns. (|40[) and 
we can solve for ifW and Wy : 

2 2 2 2 
- iA-C 2 -iA7-^C 3 r 2 + ic 3 A'- 1 exp(-2^X) (41) 



where 



W4 0) = ic 3 r 2 + C 1 T+ic 3 if- 2 exp(2XX) 



Q 



pA" 1 exp(-ATX) + 7 (42) 



, f exp(AX)|2-dX , 
*=/ ex P (A*T dX ^ 
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and p = p(u),a = a(u, X),~f = j(u) are arbitrary functions. We can now put 
constraints on the form of (141 [1 using eqn. (pt!)|) . After simpification, eqn. ([55)1 
becomes 

2C 3 K{- exp(-2KX) + 1) = (44) 
from which we conclude that C3 = (since K = is not consistent with eqn. 

69). 

Finally, we employ the remaining transformational freedom to set 7(11) = 0. 
Thus we are left with the following differential equation for H^: 

O — f) — 

2C 2 K 2 exp(2KX)a + 2C 2 K exp(2KX) 7 ^ - 2K exp(2KX) 7 ^d 

- C 2 Kcxp(2KX)a + K exp(2KX)aa - C 2 K 2 exp(2KX)T 2 

Q — 

- C\ + 2C X K 2 exp(2KX)aT + 2C x Kcxp{2K X)-^-T - C X K cxp(2KX)aT 

oX 



2CiC 2 K 2 exp(2KX)T - cxp{2KX) 



dX 



exp(KX)—?—adX - 2Kexp(2KX)a - K 2 exp(2KX)a 2 
ouoX 

da d 2 H^ d 2 H^ 
C 2 2 K 2 cxp(2KX)+exp(2KX)—a + 2-^- - 2exp(2KX)-^- 

,9ff(o) 

+ 2Kexp(2KX)— + 4K 2 exp(2KX)H {0) = (45) 

oX 

A particularly simple subcase is the case where = Wjp = = 

0, whence we obtain the special solution = Cexp(— KX), which is the 

Kaigorodov case. 

5.2 CSI example 2: 

The CSI metric can be written 

ds 2 = 2(£ 1 ® n 1 + £ 2 ® n 2 ) 

where 

t = du (46) 
n 1 = dv+[Av 2 +vH (1 \u,U,V)+H^{u,U,V)}du (47) 

+ [vVf3 + Wfj\u, U, V)]dU + [av/V + W$\u, U, V)](dV + BV 2 dU) 
t 2 = dU (48) 
n 2 = d^ + W 2 d?7, (49) 

and A, B, a and /3 are constants. 

We look for Einstein spaces in the special case: 

A = 0, a = ~2, (3 = 2B 

(where B is not specified). The metric then has the explicit form 

ds 2 = 2du®{dv+ [vH {1) +H ( - a) }du+ [2BvV + wlj ] ]dU 

+ [-2^ + W^\dV + BV 2 dU)) + 2dU ®dV + 2BV 2 dU ® df/(50) 
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The conditions for the metric to be Einstein yields the following system of 
(independent) differential equations: 

dv 2 dvou dudv 

oV z av av 

oV dv oV 

+ 8W 2 i^ - 2ff( 1 V 2 ^- - 2i?( 1 V 2 ^- + 4W 3d ^ 




du 



u 



dv du \ dv u du 



i !2V y*-2 - " u V 2 + A " - V 2 = (51) 
dudU dudV dVdU y ' 



(0) 



Y— BV 2 —V - 2BV 3 — + 2 V— 

9C/ 0V017 a^aC/ 017 

,92x^(0) Q2 W (0) dw (0) 

W v y + BV 3 Wu +B 2 V b Wv =0 (52) 
dU 2 OV 2 OV 2 [ ' 



0V 17 dU dV 

B 2 W {0) B 2 W < - 0) (9 2 iy (0) 

+ ^ 2 + w ^_^_„ s = „ (53) 

The v-dependency in eqn. (|51[) gives 
where ff« = (u,U,V), = H^(u,U,V),wij 0) = W^ 0) (u, C7, F), W^ 0) = 



14 A. Alcolado, A. MacDougall, A. Coley and S. Hervik 



Equations (|5l2j). (|53")) . (|54p do not contain any if(°) terms. The form of the 
metric ([5(1)1 is invariant under the transformations (fTCTj) and (jlip . Explicitly, 
transformation (fTU|) becomes [TO]: 



(u', v', U', V) = {u,v + h(u, U, V), U, V) (55) 

where 

H'W=H'Q-\ H'W = H'W -hHW - h u (56) 

We then define 

Wf = - 2 ± + W^-h v , (57) 

which is just rewriting the form of an arbitrary function (no u's are introduced). 
Then the other metric components transform as: 



(°) _l bt/2w(°)V - w(°) _ h _i_ m/2j, .. . nT/2fi7 (0) 



-2^+<' = (58) 



(wj; 0) + BV 2 W^) = W}j 0) - h tU + BV 2 h v + BV 2 ^ (59) 

and so the form of the metric is preserved under transformation (|10[) (with all 
arbitrary functions remaining arbitrary). 

Therefore, we use this freedom to set Wy = and obtain the solutions 

JjM = «(„) + <*> (60) 



<. =A( „,„ )t (i«i« t »a (61) 



where a, 0,^,5 are arbitrary functions. Substituting the above solutions into 
([5"Tj) gives a differential equation for 



+ 4 

V 2 dVdU 



V 2 



-ABpU + 8 , OTA 3/3 \ Jr f —ABfiU + 5 , T , 2 3/3 



2 ^7 - 77? ^ 77 + ^ 7 + 



y2 1 ' ' y3 J \ y ' 2V 2 

2 



2 y2 dU d u +2 y_L ^ =Q (62) 
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6 Conclusion 

We have presented a number of new exact 4D neutral signature VSI and CSI 
solutions that are of interest in the twistor approach to string theory [7] and 
particularly spaces admitting parallel spinors [5] . 

We note that the neutral signature case is different from the Lorentzian 
and Ricmannian cases. The Riemannian VSI and CSI cases are locally flat 
and locally homogeneous, respectively, while the Lorentzian VSI and CSI cases 
lead to the possibility of non-homogeneous Kundt spacetimes. The neutral 
signature case also leads to non- homogeneous Kundt spaces, but also allows 
for an additional "null" degree of freedom. Thus, comparing the Riemannian, 
Lorentzian and the neutral signature VSI and CSI cases, the neutral signature 
case allows for the richest variety of spaces. 
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